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Abstract. We discuss the thermal evaporation of tiny Hi clouds in the in- 
terstellar medium. Cold neutral clouds will take ~ 10 6 — 10 7 yr to evaporate if 
they are embedded in warm neutral gas and about half as long if embedded in 
hot gas. Heat flux saturation effects severely reduce the evaporation rate of cold 
neutral medium (CNM) clouds into hot gas. For CNM clouds embedded in warm 
neutral medium (WNM) the much lower conductivity results in slower evapora- 
tion. This mass loss rate could still be significant, however, if the environment 
is relatively quiescent. Partial ionization of the WNM gas would substantially 
reduce the conductivity and lengthen the lifetime of the tiny Hi clouds. The 
ultimate importance of thermal conduction to cloud evolution will depend on 
the role of turbulence and the characteristics of the medium in which the clouds 
are embedded. 



1. The Basics of Cloud Evaporation 

Cloud evaporation in the interstellar medium (ISM) is theorized to occur in re- 
gions where a cool cloud is embedded in a warmer medium. Thermal conduction 
causes heat to flow into the cloud raising its temperature and pressure at the 
boundary resulting in an evaporative outflow. The energy equation that governs 
this process is 



where A is the net cooling rate, q is the heat flux, P is the thermal pressure, v is 
the velocity and p is the density. If externally generated changes to the system 
are slow compared with the outflow speed, we may assume a steady flow. If we 
also assume spherical symmetry and subsonic outflow (v -C c, generally a good 
assumption) then equation ^ can be written as 



where m = Airpvr 2 is the (constant) mass loss rate from the cloud. As can be 
seen from this equation, cooling effectively reduces heat flux into the cloud while 
heating increases it. 

2. Thermal Conductivity in Neutral and Ionized Gas 

The heat flux can generally be expressed as q = — kVT where k is the con- 
ductivity. The conductivity depends strongly on both the temperature and 
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ionization in the gas. Generally k ~ nvl ~ T 1 / 2 /er where n is the gas density, 
u is the mean speed of gas particles, t is the mean free path and a is the col- 
lision cross section. The mean free path and collision cross section here apply 
to the carrier of the heat in the plasma - electrons in ionized gas and H atoms 
in neutral gas. For electrons in ionized gas a ~ T -2 which leads to k ~ j' 5 / 2 ) 
the classical Spitzer conductivity. In neutral gas, heat is conducted primarily by 
elas tic collisions between H atoms w ith a cross section that goes roughly as T~ - 3 
fsee iNowak <fe Ulmschneiderlll977l ) leading to k ~ T 0,8 . However, when H + is 
present, the large H + -H charge transfer cross section dramatically reduces the 
mean free path and thus the conductivity. Figure ^ demonstrates this by show- 
ing the conductivity in gas that is subject only to collisional ionization (solid 
line) compared to gas in which about a 20% H ionization fraction is maintained 
at low temperatures via an ionization rate of Th = 10 -13 s _1 (at a pressure of 
P/kB = 3000 cm -3 K) (dashed line). Clearly the ionization makes a dramatic 
difference at temperatures of T ~ 10 3 — 10 4 K such as are found in WNM/WIM 
gas. 




Temperature (K) 

Figure 1. Conductivity as a function of temperature for collisionally and 
photoionized gas. The solid line shows the conductivity for the case in which 
only collisions ionize the gas and the dashed line shows the effects of having 
a photoionization rate of 10 -13 s _1 in gas at a pressure of 3000 cm -3 K. The 
conductivity in gas in which the heat carrier is primarily neutral H is strongly 
affected by the degree of ionization in the gas because the presence of H + in 
the gas severely reduces the mean free path of the neutral H atoms. 
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3. Heat Flux Saturation 

Another important factor in limiting the heat flux under certain circumstances 
is what is known as saturation of the flux. This comes about because no matter 
how large a temperature gradient exists in the gas, the heat can only flow at the 
rate at which it can be carried by either the electrons (in ionized gas) or the H 
atoms. This limitation can be expressed as 

3 

\q\ < <7sat ~ -pc 3 , (3) 

where c is the sound speed or thermal speed of the heat carrying particles. For 
evaporation into ho t (T ~ 10 6 K) gas, sat uration is important if the saturation 
parameter, oq fsee ICowie k, McKeelll977h . is greater than 1. This parameter 
can be expressed as 

(TV/10 6 K) 3 

where R c \ is the clo ud radius and TV is the hot gas temperature (the asymptotic 
temperature in the lCowie h McKeel ( 1977ft model of steady cloud evaporation). 



Saturation, which is most important for for high TV, low P and for small clouds, 
can lead to an evaporation rate that is substantially below the "classical" rate, 

m cl (M /Myr) = 0.436(T / /10 6 K) 5 / 2 J R cl (pc). (5) 

For o"o ^ 5, which should apply to the tiny Hi clouds, the factor by which the 
mass loss rate is reduced becomes w ~ 1.22<7 5 ^ 6 . 



4. Thermal Balance and Radiative Cooling 

Radiative cooling can also act to reduce the evaporation rate for clouds. The 
strongest effects of cooling occur in the opposite regime from saturation effects 
- i.e. for high P, low TV and for large R c \. When cooling is large enough, 
evaporation is stopped altogether and the cloud is radiatively stabilized. With 
even more cooling, the cloud will condense. For a given P and T, one can 
find the critical cloud radius for stabilization - larger clouds condense, smaller 
ones evaporate. Plane parallel clouds under these assumptions should always 
condense and elongated clouds should tend to become spherical. 

Cooling via line radiation and photoelectric heating from gas and dust can 
result in a thermally bistable medium within a limited pressure range. Outside 
of that range we get only CNM (high P) or WNM (low P). This is illustrated in 
Figure |2 which makes use of data from lWolfire et all l|2003h . In this dia gram gas 
with P and n placing it at a point above the curve will tend to cool, while gas 
below the curve will heat up. Thermal conduction tends to make gas move to 
the left, possibly allowing gas to make the transition from the cold to the warm 
phase by overcoming the effects of radiative cooling. The hot phase, which would 
be off the left edge of the plot, is not a stable phase but has such a long cooling 
time, that it effectively behaves like one. The cooling lifetime of hot gas can 
be significantly reduced if its density is raised and/or its temperature reduced 
through the evaporation or mixing of colder embedded gas clouds. 
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Figure 2. Pressure vs. density for gas in thermal equilibrium in the ISM. 
The warm phase corresponds to gas along the solid portion of the curve while 
cold phase gas is for points along the dashed portion. Gas with pressure and 
density putting it on the dotted part of the curve is unstable and a small 
displacement from the curve would result in net cooling or heating, leading 
to the gas joining either the cold or warm phase. The dotted horizontal lines 
indicate the range of pressures for which a stable two phase medium could 
exist. 



Gas that is evaporated off from a CNM cloud by thermal conduction will 
have a net cooling rate that will increase sharply at first as it is heated. As its 
temperature rises it will cross the instability line in the P — n diagram and will 
have net heating. The net heating goes to zero as the temperature goes to the 
WNM temperature. If the surrounding medium is hot gas, then the temperature 
continues to increase and the net cooling increases sharply. In this case a large 
heat flux is needed to balance cooling. The net heating curve is illustrated in 
Figure 



5. Effects of the Ambient Medium on Cloud Evaporation 

CNM clouds surrounded by WNM (rather than hot gas) are protected to some 
extent from evaporation. The lower conductivity of neutral gas results in a lower 
mass loss rate. Evaporation of cold clouds into warm neutral gas is simpler than 
evaporation into hot gas in several ways. Saturation effects are not important 
and for small clouds, R c \ <^ O.f pc, cooling is also insignificant. The evaporation 
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Figure 3. Net heating rate of diffuse interstellar gas. This plot combines 
the o ptically thin cooling rate for hot gas with the net heating rate using data 
from lWolfire et all l)2003|) for cold and warm gas and shows how the cooling 
rate per unit volume sharply increases above T ~ 10 4 K. Inset is a zoom on 
the region between 100 and 10 K showing more clearly the amplitude of the 
variations in the net heating for CNM/WNM temperature gas. 

timescale is 

T „ _ SJ * x 10' <•> 

CNM clouds immersed in hot gas might be expected to evaporate quickly 
but as we discussed above, saturation effects are very important for small clouds. 
In the high o"o limit (applicable for tiny Hi clouds) we get 

— = 2.2 x 10- (J^f (^_) ( I5 ^)"' 6 yr (7) 
or expressed in terms of the cloud column 

w = « * 1° 4 (ir^r (w^f (ma) " 6 yr - (8) 

It is interesting to note that the temperature of the hot gas drops out of the 
equation for the mass loss rate in this high saturation limit. 

Depending on the surrounding medium temperature and pressure and the 
cloud characteristics, tiny CNM clouds may have short lifetimes against evapora- 
tion or relatively long ones. In Figure0]we show contours of constant evaporation 
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lifetime in the P — R c \ plane. The plot includes the case of CNM clouds em- 
bedded in hot gas (solid lines) and CNM embedded in WNM gas (dashed lines). 
Clearly the CNM can only evaporate into WNM if the two phases can co-exist, 
so the latter lines only run between the minimum and maximum pressures that 
al low for a stable two phase me dium. Also shown are lines derived from the data 
of lStanimirovic Sz Heilesl (|2005l h who give 7V(HI) and Tr for three components 
towards 3C 286. From the Figure it can be seen that if the clouds are embedded 
in hot gas, their evaporation lifetime can range from ~ 5 x 10 5 yr to ~ 2 x 10 6 
yr whereas if they are embedded in WNM the range is ~ 10 6 — 10 7 yr. 
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Figure 4. Contours of constant evaporation timescales in the R c \ — P plane. 
The dashed curves are for CNM clouds embedded in WNM (and extend only 
from P m i„ to P m ax), while the solid lines are for CNM in hot gas. The curves 
connecting the x 's are for an evaporation lifetime of 10 5 yr, while the curves 
with o's and -l-'s are for lifetimes of 10 6 and 10 7 yr respectively. The thick line 
shows the R c \—P relation for a CNM cloud embedded in WNM and radiatively 
stabilized. The dia gonal lines illustrate the R r ] — P relation for the tiny H I 
clouds observed bv iStammirovic fc Heilesl l)2005fl . The dotted lines show the 
minimum and ma ximum pressu r es for which a stable two phase medium can 
exist according to iWolfire et alJ l)2003j) . 



6. Discussion 

Recent numerical calculations show that dynamics can produce large amounts 
of gas in the "unstab l e" reg ion of the phase diagram. An example is the work of 
I Audit Hennebelld ()2005l V in which it is shown that collisions of WNM clouds 
create overpressured gas that then cools and becomes CNM gas. This is in some 
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sense the inverse process to thermal evaporation. Figure [5] shows the phase 
diagram at a time after the interaction of the two streams is well developed. 




n (cm 3 ) 

Figure 5. Pressure-density diagram color coded to indicate the gas 
mass fractions that resul t after the collision of two WNM clouds (from 
lAudit fc Hennebelld 12005 A&A, 433, 1; reprinted by the kind permission 
of the authors and Astronomy & Astrophysics. See original article for full 
color version of this figure). The parallel diagonal lines are isothermal curves 
for T = 5000 K and T = 200 K. The full black line is the therm al equilibrium 
curve. (Note that this curve is somewhat different from the I Wolfire et alJ 
(2003) equilibrium curve). Much of the mass has been compressed so as to 
be out of equilibrium and subject to cooling. In this way gas is fed from the 
WNM to the CNM as a result of the cloud dynamics. 



Such results lead to the question of the relative importance of thermal 
conduction/evaporation versus dynamical processes such as this sort of triggered 
condensation. It is clear that the evaporation timescale is long enough that 
evaporation is unlikely to be important in regions experiencing such strongly 
dynamical processes. On the other hand in more quiescent regions, evaporation 
could dominate. The overall importance of conduction thus depends on an 
understanding of the larger scale evolution of the ISM and such things as the 
frequency of shock crossing for a given shock strength or the frequency and 
relative speed for cloud-cloud collisions. The old questions of the filling factor of 
CNM, WNM, WIM (warm ionized medium) and HIM (hot ionized medium) and 
their relative arrangement in the ISM clearly also has bearing on the lifetime 
for tiny H I clouds. As we have shown, the ambient medium in which the clouds 
are embedded could be crucial to determining their lifetimes. 
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7. Conclusions 

Though many may have the impression that small clouds, such as the tiny 
Hi clouds, ought to evaporate very quickly in hot gas, the effects of heat flux 
saturation slow the process substantially. Even without appeal to tangled mag- 
netic fields the lifetime against evaporation of these clouds should be on the 
or der of 10 6 yr. If th e clouds are instead surrounded by warm neutral gas, a 
la iMcKee fc Ostrikerl (Il977ft . their lifetimes against evaporation should be in- 



creased, though only to <J 10 7 yr. Thermal evaporation should occur within 
neutral gas, i.e. between CNM and WNM, and is simpler than between CNM 
and hot gas since neither the magnetic field nor radiative cooling nor saturation 
should play much of a role. The one complication in this case is that partial 
ionization of the warm gas can substantially reduce its thermal conductivity. 
More study is needed in this area, particularly regarding possible observational 
diagnostics of such warm-cold evaporation. The final word on the importance 
of thermal conduction to the evolution of the tiny H I clouds will have to await 
a more thorough understanding of the environments in which they exist and the 
overall morphology and dynamics of the diffuse ISM. 
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